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Abstract 

The atmosphere of Mars significantly attenuates the heavy ion component of the 
primary galactic cosmic rays (GCR), however increases the fluence of secondary light ions 
(neutrons, and hydrogen and helium isotopes) because of particle production processes. 
We describe results of the quantum multiple scattering fragmentation (QMSFRG) model 
for the production of light nuclei through the distinct mechanisms of nuclear abrasion and 
ablation, coalescence, and cluster knockout. The QMSFRG model is shown to be in 
excellent agreement with available experimental data for nuclear fragmentation cross 
sections. We use the QMSFRG model and the space radiation transport code, HZETRN to 
make predictions of the light particle environment on the Martian surface at solar 
minimum and maximum. The radiation assessment detector (RAD) experiment will be 
launched in 2009 as part of the Mars Science Laboratory (MSL). We make predictions of 
the expected results for time dependent count-rates to be observed by RAD experiment. 
Finally, we consider sensitivity assessments of the impact of the Martian atmospheric 
composition on particle fluxes at the surface. 



INTRODUCTION 


Human beings are scheduled to visit Mars in the 2030 to 2040 time-period. Prior to their 
arrival an accurate assessment of the radiation environment on the Mars surface is 
required. Particle flux distributions, which describe particle type (charge and mass), j, 
kinetic energy, C(MeV/u) and particle directions are necessary for a complete description 
of radiation fields in space because the use of dose equivalent carries large uncertainties in 
projecting risk of late effects from galactic cosmic rays (GCR) [1-3]. The primary GCR 
are well described by empirical models based on extensive measurements [4]. The 
description of the GCR on the planetary surface of Mars is expected is more challenging 
because of larger amount of material ions traverse in the Mars atmosphere and backscatter 
of neutrons and other light ions from the Mars regolith [5], Nuclear fragmentation and 
energy loss processes are the main physical mechanisms in GCR transport. In this report 
we briefly summarize a quantum multiple scattering model of nuclear fragmentation 
(QMSFRG) and describe its extension to incorporate nuclear coalescence as a mechanism 
of light ion production. The RAD experiment will be launched in 2009 as part of the Mars 
Surface Lander (MSL) [6]. In this report we discuss predictions of the expected results for 
time dependent count-rates to be observed by RAD experiment on Mars, and evaluate the 
role of the altitude of a landing site and the atmospheric composition. 

METHODS 

Quantum Fragmentation Model 

The QMSFRG considers the multiple scattering series for two heavy ions and introduces 
the impulse and eikonal approximation for the total momentum transfer vector in order to 
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obtain a closed-form solution to the abrasion cross section excitation spectrum [7-10]. For 
inclusive reactions where a single fragment originating in the projectile is measured, 
closure is perfonned on the final target state with a momentum vector denoted p x used to 
represent these states. The reaction is assumed to proceed through the abrasion step 
producing a pre-fragment, F* and fireball piece from the projectile-target overlap denoted, 
R, followed by the ablation step where the final projectile fragment, F is formed after 
nuclear de-excitation. The total momentum transfer is q = p r - p x where p T is the initial 
target momentum. The pre-fragment, F * excitation spectrum following nucleon or alpha 
particle abrasion was found [7,8] as an impact parameter dependent convolution of the 
pre-fragment excitation response for a transition of the pre-fragment core from state n to 
n' and the project fireball response, average over the target ground-state wave function, 
I T>, 
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where b (b ) is the impact parameter, and q the momentum transfer. The abrasion response 
is defined as the quantum matrix elements for the interaction of the projectile fireball, R 
with the target, T after performing closure over the final fireball states (denoted | R>) 
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where the Qrt represent the fireball-target profile operator, and kR the projectile fireball 
momentum vector. The abrasion-response represents a complicated many-body operator 
that is solved by approximation using closure over the target and fireball states for 
evaluating the pre-fragment distribution. The one-particle abrasion response has been 
evaluated using the shell model response functions [8, 9]. The pre-fragment excitation is 
described in terms of the transition matrix 

(3) P n n ,(b,b' ,E pt ) = <F\\ Q F A(b') \ F >< F | Q F . T (b) \ F *„ > 
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where quantum matrix elements for the pre-fragment excitation are evaluated over the 
many-body profile operators, Qf*t- In the model a convolution approach is used to derive 
the multi-nucleon abrasion excitation spectrum from the single-fragmentation term [7,8]. 
The resulting excitation spectrum is broad with a shape similar to a log-nonnal distribution 
with mean energies from 20 to 30 MeV for one-nucleon removal. 

The de-excitation of the pre-fragments in nuclear ablation is described in a stochastic 
process using a Master equation for de-excitation by particle emission [9,10], If fb(E,t) is 
the probability of finding the nuclei b at time t with excitation energy Eb and Pb,k(E) be the 
probability that the nuclei, b will emit ion k with energy E, then the Master equation is 

(4) dUE l f ) = Y.\ dE fM^PJ.E)-^\dEf b (El,bP h . k (E) 

dt ] k 

The first-term on the right side of Eq. (4) corresponds to gains by decays of ion a emitting 
ion j to form ion, b, and the second tenn from losses due to decays of ion b into ion c by 
emitting ion k where the j (or k) are light-particle emissions (n, p, d, t, h, or a). The 
probability of finding the nuclei b at time t with E\ can be divided into stable and unstable 
parts depending on the lowest excitation energy of E * denoted min[ S bj ], 

(5) f(El,t) = g(E* b ,t) + h(E* b ,t) 

where g and h denote the stable and unstable parts of f respectively. As t — >co, we have 
h — >0 such that 

lim * * 

(6) f(E b ,t) = g(E b ,co) 
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The probabilities for a single-step decay of nucleus a, are defined as 
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with a stable daughter nuclei formed corresponding to the probability 


< 8 > - I 


E -S„ 




P,(E)dE 


and an unstable daughter nuclei formed with probability 
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These probabilities obey 

(10) F aJ = G aj + H aj 

and the normalization condition, '^ j F a . = 1 . Because only stable nuclei remain at 

j 

long-times, the differo-integral Eq (4) is separated using the above definitions as 

(11) «,(£>) = 


and similarly 
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In evaluating fragmentation cross sections, equations (11) and (12) are solved by iteration 
up to excitation energies of 200 MeV [9,10], and using the approximation of Campi and 
Hufner [11] at higher values. The use of the statistical decay model with accurate nuclear 
level densities that include nuclear shell structure effects and the use of measured values 
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for the nuclear masses are important for accurately predicting the odd-even effects in the 
fragment spectrum. 

To estimate nuclear coalescence contributions to d, t, h, and a production we consider the 
momentum distribution for proton (p ) and neutron (n) production in the QMSFRG model 
[10, 13], and apply the usual coalescence model [12], which consists of an A n -fold folding 
of the p and n- distribution to form a light ion of mass number A n . The n and p production 
cross sections are then reduced by the appropriate balancing of the cross sections resulting 
in light nuclei coalescence. In these estimates only the contributions from abrasion 
momentum distribution’s is used. 

Space Radiation Transport Code 

We use the HZETRN code of Wilson et al. [14, 15] in our predictions of the Mars surface. 
HZETRN solves for the spectrum of heavy ion fragment from projectile and target nuclei 
in the continuous slowing down and straight-ahead approximations. Integration over the 
Martian atmospheres is perfonned as described previously [3], We have modified the 
HZETRN code in several ways: First we have extended the ion grid to 190-ions from the 
earlier 59- [13] or 170-ion [16] grids with the ions selected after studying the fragment 
spectra for most GCR nuclei for primaries of Z=1 to 28. Second we use the QMSFRG data 
base as described above. Finally, we use the primary GCR model of Badhwar et al. [4], 
but re-distribute the most abundant element used in this work amongst the isotopic 
composition found for the primary GCR as described in Cucinotta et al. [16]. 

RESULTS AND DISCUSSION 

We first describe comparison of the QMSFRG model to heavy ion fragmentation cross 
sections. Our recent publication [16] made extensive comparisons of QMSFRG to 
experimental cross sections for Ne, Mg, Ar, and Fe. In Figure 1 we show comparisons of 
the correlations between QMSFRG and the recent data by Zeitlin et al. [17] for “ Si 
fragmentation cross sections at several energies for C and Al targets. As in our earlier 
work, QMSFRG demonstrates good agreement with fragmentation cross sections for 
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projectile fragments as function of target mass and beam energy with over 85% of the 
cross sections agreeing to within +25 percent of measurements. Light ions are produced 
by several mechanisms including nuclear abrasion, ablation, coalescence, cluster 
knockout, and electro-magnetic dissociation. In Table 1 we show contributions for several 
of these mechanisms for 16 0 fragmentation on several targets and compare to the 
experimental data of Olsen et al. [18]. Coalescence is predicted to be a dominant 
mechanism for d, t, h, and a-particle production. Calculations include cluster knockout 
only for a-particle production and appear to be most important for H and other light 
targets. Studies of d, t, and h cluster knockout maybe warranted, especially for target 
fragmentation data bases. 

In Table 2 and Figures 2 and 3 we show predictions for the Mars surface. Table 2 shows 
the elemental fluence at solar minimum where the solar modulation parameter, ® of 428 
MV is used [4] and for the year 2010 a modulation parameter of 1000 MV. The Mars 
atmosphere significantly depletes the HZE components and the Z=l, 2 ions will be the 
main inducers of biological risk for astronauts occupying the Mars surface. The angular 
view of RAD will be no more than 90 degrees because of the mass limitations on MSL, 
and we show results for the complete horizon and a restricted vertical horizon of 60 and 90 
degrees in Table 2. In Figure 2 we show results for the energy spectra and average 
fluence of light and heavy ions for different altitudes relative to the mean surface. The 
shape of the spectrum below 200 MeV/u changes dramatically due to production from 
GCR heavy ions or target atoms in the atmosphere. We have not included the backward 
neutron components in these comparisons, which significantly increase the neutron 
components below about 100 MeV and can lead to increases on more than one order of 
magnitude for neutrons below about 10 MeV especially for landing sites on Mars regolith. 

The Mars atmosphere is largely CO 2 , however contains a small amount of other 
constituents [19-21]. It is also of interest for improved understanding of the possible 
history of life on Mars to consider the evolution and possible changes of the Martian 
atmosphere and their impact on radiation effects for possible life in the past. Figure 3 
makes predictions for several possible atmospheric conditions including a methane 
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atmosphere or mixed CH 2 and CH 4 atmosphere. Along with the vertical height of the 
atmosphere, the mass contribution of hydrogen makes the most significant modulation of 
the surface radiation environment. The relative contributions from non-hydrogen atoms 
has only a small effect on the charged-particles observed on the Mars surface, however are 
expected to have a slightly larger dependence for the albedo neutrons [5]. 

In summary, we have shown that the quantum multiple scattering model of nuclear 
fragmentation (QMSFRG) provides an accurate data base for GCR transport problem 
applications. Because the atmosphere of Mars represent one to two mean free paths for 
GCR heavy ions leading to buildup of light ion fragments, we have made predictions of 
the flux rates to be expected to be observed by RAD in 2010 for light and heavy particles, 
which include the nuclear coalescence contribution to the cross sections. Because the 
integral flux of GCR nuclei changes slowly with atmospheric parameters or 
geographical location of Mars, measurements of charge-specific energy spectra 
between 10 and 200 MeV/u are shown to provide the optimal tests of radiation 
transport codes. Future work will aim to improve the calculation of cross sections in the 
coalescence and cluster knockout models [8,12], include the albedo neutron and light ion 
contributions, and describe detailed geometric models of the RAD configuration on MSL 
in preparation for the data analysis phase of the experiment. 
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FIGURE CAPTIONS 


Figure 1. Comparisons of fragmentation cross sections in QMSFRG model to 

28 

experimental data [19] for Si fragmentation on C and A1 targets at several energies. 


Figure 2. Upper panel shows calculations of the energy spectra of Z=1 and 2 ions incident 
on Mars and the average mean altitude projected for 2010 (<f>= 1 000 MV). The lower panel 
shows calculations of the total annual flux of Z= 0, 1, 2, 6, and 8 ions for different mean 
altitudes on Mars projected for 2010 (0=1000 MV). 


Figure 3. Calculations of the dependence of the Mars surface point dose equivalent on 
atmospheric composition and density. Different atmospheric mixtures are shown including 
the temptative detection of methane (CH4), which vary between lOppb (parts per billion) 
and a maximum of 300ppb of CH4. Two theoretical end cases of a pure C0 2 and pure CH 4 
atmosphere are shown for reference.. 
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Figure 1. 
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Table 1A. Contributions to light ions production in I6 0 fragmentation on H at 2.1 GeV/u in QMSFRG compared to experiment [19]. 
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Table 2a. Projections of Elemental GCR fluences for RAD Surface operations in 2010 
(0=1000 MV). 


Charge Group, Z 

Fluence, 1/(cm 2 hr) 
All angles 

Fluence, 1/(cm 2 hr) 
90 degree viewing 
cone 

Fluence, 1/(cm 2 hr) 
60 degree viewing 
cone 

0 

13841 

4374 

2751 

1 

11576 

5756 

3817 

2 

425 

288 

198 

3 

3.32 

2.30 

1.58 

4 

2.40 

1.66 

1.14 

5 

3.23 

2.38 

1.66 

6 

7.08 

5.40 

3.82 

7 

2.34 

1.77 

1.25 

8 

4.46 

3.56 

2.56 

9 

0.30 

0.23 

0.16 

10 

0.80 

0.64 

0.46 

11 

0.31 

0.24 

0.17 

12 

0.78 

0.64 

0.46 

13 

0.20 

0.16 

0.12 

14-17 

0.42 

0.35 

0.26 

18-23 

0.23 

0.19 

0.14 

>23 

0.18 

0.16 

0.12 


Table 2b. Same as Fig 2a near Solar Minimum (0=428 MV) 


Charge Group, Z 

Fluence, 1/(cm 2 hr) 
All angles 

Fluence, 1/(cm 2 hr) 
90 degree viewing 
cone 

Fluence, 1/(cm 2 hr) 
60 degree viewing 
cone 

0 

24779 

8131 

5127 

1 

22552 

11711 

780 

2 

770 

527 

364 

3 

5.94 

4.19 

2.89 

4 

4.20 

2.96 

2.04 

5 

5.58 

4.17 

2.93 

6 

12.09 

9.37 

6.67 

7 

3.89 

2.99 

2.12 

8 

7.40 

5.99 

4.34 

9 

0.48 

0.36 

0.26 

10 

1.25 

1.01 

0.73 

11 

0.48 

0.38 

0.27 

12 

1.19 

0.99 

0.72 

13 

0.31 

0.25 

0.18 

14-17 

0.96 

0.80 

0.59 

18-23 

0.33 

0.28 

0.20 

>23 

0.25 

0.22 

0.17 
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